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ABSTRACT 

We report the results of mapping observations of the bipolar nebula with 
SiO maser emission, IRAS 19312+1950, in the CO (J = 1-0 and J = 2-1), 
13C0 (J = 1-0 and J = 2-1), C^^O (J = 1-0), CS (J = 2-1), SO (Jk = 83- 
2i) and HCO"'" (J = 3-2) lines with the Berkeley-IUinois-Maryland Association 
array. Evolutional status of this source has been evoking a controversy since 
its discovery, though SiO maser sources are usually identified as late-type stars 
with active mass loss. In line profiles, two kinematical components are found: a 
broad pedestal component and a narrow component. Spatio- kinetic properties 
of a broad component region are roughly explained by a simple spherical outflow 
model with a typical expanding velocity of an AGB star, though some properties 
still conflict with those of a typical AGB spherical outflow. A narrow component 
region apparently exhibits a bipolar flow. The angular size of the narrow compo- 
nent region is spatially larger than that of a broad component region. Intensity 
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distribution of the CS emission avoids the central region of the source, and that of 
an SO broad component emission exhibits a small feature peaked exactly at the 
mapping center. According to the present results, if a broad component really 
originates in a spherical outflow, an oxygen-rich evolved stellar object seems to 
be a natural interpretation for the central star of IRAS 19312+1950. 

Subject headings: masers — circumstellar matter — stars: imaging — stars: 
individual (IRAS 19312+1950) — stars: winds, outflows — ISM: jets and out- 
flowsmasers — circumstellar matter — stars: imaging — stars: individual (IRAS 
19312+1950) — stars: winds, outflows — ISM: jets and outflows 

1. Introduction 

Evolutional status of the SiO maser source, IRAS 19312+1950, has been evoking a 
controversy since its discovery in our SiO maser survey toward a sample of relatively cold 
IRAS sources (Nakashima Sz Deguchi 2000, 2003). Infrared colors and morphology of the 
source are fairly reminiscent of those of asymptotic giant branch (AGB) stars (or post-AGB 
stars); in fact, SiO maser sources are usually identified as oxygen-rich (0-rich) latc-type 
stars with active mass loss (see, e.g., Habing 1996). Some properties of the source, however, 
arc consistent with those of young stellar objects (YSOs) in molecular clouds (Nakashima 
& Deguchi 2000, 2004; Nakashima et al. 2004; Deguchi et al. 2004). IRAS 19312+1950 
exhibits a twin-peaked SiO maser fine profile (Nakashima & Deguchi 2000); this exhibits 
strong time variation (Deguchi et al. 2004). The twin-peaked SiO maser profile is rarely 
seen in AGB stars, and is reminiscent of the SiO maser profile of the Orion KL star-forming 
region (Snyder & Buhl 1974; Moran et al. 1977; Plambcck et al. 1990; Doeleman et al. 1999). 
H2O (22.2351 GHz) and OH (1612, 1665, and 1667 MHz) maser lines have also been detected 
(Nakashima & Deguchi 2000, Lewis 2000 in private communications). The H2O maser line 
profile shows several intensity peaks with strong time variation (Deguchi et al. 2004). The 
OH maser profiles show complicated spiky profiles, which are unusual in AGB stars, but the 
line intensity ratios between three different OH maser transitions are consistent with those 
of type II OH maser sources in AGB stars (i.e., the 1612 MHz sateUite line is the strongest 
among the three fines, see, e.g., Habing 1996). 

A rich set of molecular species, detected toward IRAS 19312+1950, has recently raised 
an exotic possibility that this source is an AGB/post-AGB star lying in a dark cloud (Deguchi 
et al. 2004), though the idea is not conclusive yet. In fact, 16 molecular species (including 
isotope species) have been detected so far toward the source (Deguchi et al. 2004); this rich- 
ness of molecular species is reminiscent of that of OH 231.8+4.2 (rotten egg nebula, see. 
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e.g., Morris et al. 1987). The molecular line profiles often show two kinematic components: 
strong lines with small line-widths of ~ 2 km s~^ (narrow component) and weak-symmetric 
lines with large line-widths of > 60 km s~^ (broad pedestal component) (Nakashima et al. 
2004; Nakashima & Deguchi 2004). The broad pedestal component seems to be attributed to 
an AGB wind in terms of its large line-width, but the narrow component is rather difficult 
to explain simply by an AGB/post-AGB star origin (Nakashima & Deguchi 2004). Such 
a narrow line profile is relatively rare but known in several late-type stars (for example, 
EP Aqr, X Her, EU And, and EM Gem; Knapp et al. 1998; Kerschbaum & Olofsson 1999; 
Jura & Kahane 1999; Winters et al. 2003; Nakashima 2005). A recent interferometric radio 
observation of IRAS 19312+1950 in the HCO+ J = 1-0 fine (Nakashima & Deguchi 2004) 
has revealed existence of a bipolar flow and a (conceivably) spherical outflow with a small 
angular size. If IRAS 19312+1950 is a true AGB/post-AGB star, as suggested in our pre- 
vious papers, the rich molecular environment of the source offers a compelling opportunity 
to investigate the chemical environment of a unique 0-rich AGB/post-AGB envelope. Oth- 
erwise, even if the source identifled as a YSO, it must provide an alternative opportunity to 
investigate a YSO with SiO masers, because the number of YSOs, emitting SiO masers, are 
very limited. However, further astrophysical interpretation of the source is contingent upon 
secure identiflcation of its evolutional status. 

In this paper, we present the results of a set of high resolution radio interferometric 
observations in 8 molecular rotational fines with the Berkeley-IUinois-Maryland Association 
(BIMA) array. Main objectives in the present research are (1) to reveal a precise spatio- 
kinetic properties of narrow and broad component regions, and (2) to investigate chemical 
properties of the source, with high-angular-resolution radio mapping observations. We make 
use of the present results to consider the evolutional status of the source. The outline of 
this paper is as follows: in section 2, details of observations and results are presented. In 
section 3, spatio-kinctic properties of narrow and broad component regions arc discussed; 
chemical properties of the source are also briefly discussed in this section. Finally, results 
are summarized in section 4. 



2. Observations and Results 

Interferometric observations of IRAS 19312+1950 were made with the BIMA array; 
the instrument has been described in detail by Welch et al. (1996). We observed the GO 
(J = 1-0 and J = 2-1), I^CO (J = 1-0 and J = 2-1), G^^O (J = 1-0), GS (J = 2-1), SO 
( Jk = 32-2i) and HCO+ (J = 3-2) fines in two observing periods: February - June 2003 (for 
CO) and December 2003 - January 2004 (for CS, SO and HCO+). Rest frequencies of the 
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lines were taken from Lovas (1992). The observations were interleaved every 25 minutes with 
the nearby point source, 1925+211, to track the phase variations over time. The absolute 
flux calibration was determined from observations of Uranus and MWC349, and is accurate 
to within 20%. Typical single sideband system temperatures ranged from 200 K to 300 K 
(at A ~3 mm) and from 700 K to 1500 K (A ~1 mm). The velocity coverage was about 380 
km s~^ (A ~3 mm) and 150 km s~^ (A ~1 mm), using three different correlator windows 
with a bandwidth of 50 MHz each. The velocity resolution was about 1 km s~^ (A ~3 mm) 
and about 0.5 km s^^ (A ~1 mm). The phase center of the maps is 19'^33'"24.4^ 19°56'54.8" 
(J2000) corresponding to the IRAS position of the target. Data reduction was performed 
with the MIRIAD software package (Sault et al. 1995). Standard data reduction, cahbration, 
imaging, and deconvolution procedures were followed. The rms noises, integration times on 
the source, beam sizes given by robust weighting of the visibility data, position angles of 
the beams, and array configurations are summarized in Table 1. Note that about 50% of 
the visibility data of the CO J = 1-0 line observations were fiagged out because of severe 
mechanical problems on a correlation system; therefore, the integration time of the CO 
J = 1-0 line observations is practically a half of the time indicated in Table 1. 



2.1. Spectra 

In Figures 1 we present spectra of molecular lines observed, averaged spatially over a 
circle with an angular diameter of 15"; this averaging circle is centered on the mapping center. 
The line profiles, shown in Figure 1, are rather similar with those taken in previous single- 
dish observations with careful background subtraction (Nakashima et al. 2004; Deguchi et 
al. 2004); this suggests that almost all of the single-dish flux is recovered by the present 
interferometric observations. In comparison with the previous single-dish data taken by 
the Nobcyama 45-m telescope (Nakashima et al. 2004), about 70% of flux, detected by the 
Nobeyama 45-m telescope, is recovered in the present interferometric observations; the main 
reason of the difference of flux is due presumably to difference of calibration methods in each 
telescope. The line profiles exhibit two kinematical component, as reported by Nakashima et 
al. (2004): "narrow component" peaked strongly at Visr =38-40 km and "broad pedestal 
component" . The narrow component is seen in all lines observed, but the broad component 
appears to be limited in the CO J = 1-0, CO J = 2-1, ^^CO J = 2-1, and SO Jk = 32-2i 
lines (the broad component is particularly prominent in the CO J = 1-0 and CO J = 2-1 
lines). The ^^CO J = 1-0 and CS J = 2-1 fines seem to have a weak shoulder at the 
blueshifted side of the narrow component. According to intensity maps, the shoulder feature 
of the ^^CO J = 1-0 line seems to be a part of the broad pedestal component, whereas 
that of CS J = 2-1 seems to be a part of the narrow component (i.e., narrow and broad 
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components show different spatio- kinetic properties; we discuss this in later sections). Our 
previous observation in the HCO"*" J = 1-0 hne shows a broad pedestal component with 
a parabolic profile (Nakashima & Deguchi 2004); however, the present observation in the 
HCO"*" J = 3-2 very marginally shows a weak broad component. This is due presumably 
to high noise level at the frequency very close to an edge of band coverage. Peak velocities 
of narrow components of the CO J = 1-0 and CO J = 2-1 lines are lying at ~ 40 km 
s~^, while narrow components of other lines are peaked at ~38 km s~^. This disagreement 
of peak velocities is due to strong absorption features lying at Visr ~35 and ~45 km s~^. 
A single-dish spectrum in the CO J = 1-0 with no background subtraction (see, Fig. 1 
in Nakashima et al. 2004) shows a strong emission feature at V\si ~35 km s~^; therefore, 
the absorption feature at V\si ~35 km s~^ is most likely real absorption due to foreground 
material, detected by the single-dish observation as an emission line. Although no strong 
emission was detected at V^sr ~45 km in the single-dish observation, the 45 km 
feature also seems to be attributed to foreground clouds. This is because, as we will see 
in channel- velocity maps later, red- and blueshifted sides of the absorption feature lying at 
Visr ~45 km exhibit almost same spatio-kinematic properties. The ^^CO J = 1-0 and 
C^^O J = 1-0 lines, as well as the ^^CO lines, seem to exhibit weak absorption features. 
Nakashima et al. (2004) suggested, on the basis of their single-dish data, that a full width of 
the broad pedestal component of the CO J = 1-0 line reaches, at least, to 60 km s~^; in fact, 
the present intcrfcrometric data show that the full-width at the zero intensity level reaches 
nearly to 80 km s-\ The CO J = 2-1, ^^CO J = 2-1 and HCO+ J = 3-2 lines were first 
observed toward IRAS 19312-1-1950 in the present observations. Intensity of the CO J = 2-1 
and ^^CO J = 2-1 lines is clearly larger than that of the CO J = 1-0 and ^^CO J = 1-0 
lines. Approximate intensity ratios between higher and lower transitions of CO lines are [CO 
(J = 2-l)]:[C0 (J = 1-0)] = 1.76:1 (narrow component), [CO (J = 2-l)]:[C0 (J = 1-0)] = 
8.5:1 (broad component), [^^CO (J = 2-l)]:[i3co (J = 1-0)] = 6.7:1 (narrow component). 
The intensity ratios of narrow components were calculated with peak intensity. The ratios 
of broad components were calculated with integrated intensity, which was integrated over 
particular velocity ranges (i.e., 15-20 km s~^). The ratio of a broad component of the ^^CO 
lines is not given here, because that is too weak to measure fluxes in reliable accuracy. In 
the C^^O J = 1-0 line proflle, a weak emission-like feature is seen at Visr ~ HO km s^^. We 
verified, with integrated intensity maps, that the possible high velocity component at V^gr ~ 
110 km s~^ is not real. 

Continuum emission was not detected at a 3 cr rms level both in 1 mm and 3 mm bands; 
to confirm the continuum emission, frequency channels were integrated over 150 MHz using 
the opposite-side bands to line observations. Upper limits of the continuum radiation flux 
at 3 mm and 1 mm bands (center frequencies of integration ranges are ~107 GHz and ~228 
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GHz) are 2.2 x 10~^ Jy and 4.9 xlO~^ Jy, respectively. 

2.2. Intensity Maps 

Figures 2 and 3 show velocity channel maps of the ^^CO J = 1-0 and ^^CO J = 2-1 lines, 
respectively. The structure seen in Figures 2 and 3 is clearly resolved by our synthesized 
beam. Broad and narrow components mentioned in section 2.1 exhibit clearly different 
spatial structure in the velocity channel maps; this difference is prominent especially in 
Figure 3. In Figure 3, the velocity channels of 35-37 km s~^ roughly correspond to a 
narrow component, and those of 24-32 and 40-46 km s~^ correspond to a broad component. 
Broad component features appear to be roughly spherical morphology, and the sizes of the 
broad component features are smaller than those of narrow component features. The narrow 
component features extends in the north-west to south-east direction. In Figure 2, a narrow 
component is predominantly seen in velocity channels of 29-30 km s~^ and 35-38 km s~^. No 
strong features are seen in channels of 32-34 km s~^ due to absorption in the corresponding 
velocity range. Broad component features are somewhat difficult to see in Figure 2 because 
of its weak intensity; however, a small feature peaked at the mapping center is seen in 
the 40 km s~^ channel. In both of Figures 2 and 3, narrow component features appear 
to exhibit systematic variation as a function of radial velocity: narrow component features 
symmetrically vary with respect to the velocity ~36.5 km s^^. In maps of a blueshifted side 
(i.e., 35 and 36 km s~^ channels), narrow component features extend mainly to the north-west 
direction, and exhibit almost no extension to the south-east direction; in contrast, features 
at a redshifted side (37 and 38 km s~^ channels) extend mainly to the south-east direction, 
and almost no extension is seen to the north-west direction. Similar variation of structure is 
seen also in velocity channel maps of the C^^O emission; the maps are presented in Figure 
4. 

In Figures 5 and 6 we present velocity selected channel maps of the CO J — 1-0 and 
CO J = 2-1 lines, respectively. As well as the ^^CO and C^^O hues, spherical and elongated 
features are seen: velocity channels of 37-39 km s~^, corresponding to a narrow component, 
exhibit largely extended features elongated in the north-west to south-east direction, and 
those of Visr < 26 km s~^ and V^sr > 42 km s~^, corresponding to a broad pedestal component, 
exhibit small sizes and apparent sphericity. A narrow component region of the CO J = 1- 
emission is about 2 times larger than that of the CO J = 2-1 emission; in contrast, a 
broad component region of the CO J = 1-0 emission appears to be about 2 times smaller 
than that of the CO J = 2-1 emission, although signal-to-noise ratios in the CO J = 1-0 
observations are somewhat small because of significant elimination of bad data as mentioned 
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in the previous section. The systematic change of structure with respect to Visr ~36.5 km 
s~^, seen in the ^^CO and C^^O hnes, is not seen in the CO J = 1-0 and CO J = 2-1 hues 
due presumably to interruption of absorption features. 

To see overall intensity distribution of broad and narrow components, in Figures 7 and 
8 we present velocity integrated intensify maps of each kinematic component. Figure 7 

shows integrated intensity maps of the narrow components. To clarify a systematic change 
of structure with respect to l^sr ~ 36.5 km s^^, we integrated intensity of narrow components 
over two different velocity ranges: 35-36 and 37-38 km s~^. In Figure 7 intensity maps of 
these integration ranges arc indicated by the thick (35-36 km s~"^) and thin (37-38 km s~"^) 
contours. The thick and thin contours exhibit clear difference in their maps; this difference 
is prominent in the ^^CO J = 1-0, ^^CO J = 2-1 and C^^O J = 1-0 lines. At a glance, 
features composed by thin and thick contours are reminiscent of a bipolar flow. In the case 
of ^^CO, a feature extends in the north-west to south-east direction, and has an angular size 
of about 50" X 30" with a position angle of about 130°. These angular sizes correspond to 
linear sizes of 1.5 x 10^* and 9.0 x 10^'' cm at 2 kpc; similarly, 3.0 x 10^* and 1.8 x 10^^ cm 
at 4 kpc (see, e.g., Nakashima et al. 2004, on distance estimation to the source). Position 
angles of major axes of the elongated features seem to be slightly varies from line to line. 
If this feature is of a bipolar out-flow, the north-west side should be closer to us, and the 
south-east side should be further from us. In the cases of ^^CO J = 1-0 and J = 2-1 
lines, emission in the range of Visr = 35-36 km is almost totally obscured by absorption 
features; therefore, features are almost absent in thick contours. Although features in Figure 
7 are mainly elongated in the north-west to south-east direction, another minor elongation 
is seen in the south-west side of the mapping center especially in maps of the ^^CO J = 1-0 
and C^^O J = 1-0 lines; near-infrared structure (sec, e.g. Nakashima & Deguchi 2004) is 
also elongated in this direction. Note that intensity peak positions of the ^^CO and C^^O 
narrow components are lying at 3" north-west of the center, and do not correspond to the 
mapping center,. 

In Figure 8 we present velocity integrated intensity maps of broad components. The 
right and left columns represent maps of blue- and redshifted wings of broad components, 
respectively. In addition to CO maps, in Figure 8 a similar map of the SO line is also 
presented. Features seen in Figure 8 are clearly resolved by our synthesized beam except for 
SO, and morphology of features are totally different from that of narrow components. At a 
glance, structure seen in Figure 8 seems to be roughly spherical. Intensity peak positions of 
red- and blueshifted wings mutually correspond within a beam size. The angular diameter 
of the ^^CO (largest one) J = 2-1 spherical feature is about 11", corresponding to 3.3 x 10^^ 
at 2 kpc (or 6.6 x 10^'' cm at 4 kpc). In the case of SO, the structure is not resolved by our 
synthesized beam; this suggests that angular sizes of the SO broad component regions is less 
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than about 3". 

In Figure 9, we present velocity integrated intensity maps of narrow components of CS, 
SO and HCO+. A map of the HCO+ J — 1-0 line is also presented for comparison; the 
HCO+ J = 1-0 data were taken from Nakashima & Deguchi (2004). The intensity maps are 
superimposed on a near-infrared J, X-bands composite image. In Figure 9, similarly 
with Figure 7, the thick and thin contours represent velocity integration ranges of 34-36 
and 37-39 km s~^. Note that scales of the maps are different from those of CO maps (i.e., 
the region of Figure 9 are about 6 times smaller than that of Figures 2-5). In contrast to 
the CO lines, no significant structure variation with respect to radial velocity is seen in the 
CS, SO and HCO+ lines, and structure of these lines is clearly different from that of the 
CO line. Features seen in Figure 9 appear to avoid the mapping center, where near-infrared 
emission exhibits a bright peak. A notable intensity peak is lying at about 3" north-west of 
the mapping center; this position is common in all four lines, and corresponds with intensity 
peak positions of the ^^CO and C^^O narrow components. As mentioned in the section 2.1, 
the CS line profile shows a weak shoulder at the redshifted side of its peak; this shoulder 
component shows an intensity peak at the position 3" north-west of the mapping center. The 
SO map exhibits another weak peak at about 2" south-east of the mapping center; HCO"'" 
(J = 3-2) exhibits a peak at about 7" east of the mapping center. A notable point in Figure 
9 is asymmetry of intensity distributions: all lines appears to be stronger at the north-west 
side than at the south-east side of the mapping center. This asymmetry seems to show 
significant difference of physical condition between the north-west and south-east sides. 

3. Discussion 

In this section, we have discussed evolutional status of IRAS 19312-t-1905 and related 
problems by making use of the present results. 

3.1. Properties of the Broad Component Region 

A prompt interpretation of a broad component region, apparently exhibiting spherical 
morphology in intensity maps, would be a spherically expanding outfiow. For better under- 
standing of spatio-kinetic properties of the broad component region, in Figure 10 we present 
a structure size vs. radial velocity diagram using the CO J = 1-0, CO J = 2-1 and ^^CO 
J — 2-1 data. Though the SO line also exhibits a weak broad component, here we do not 
use the SO data, because structure of the SO broad component region is not resolved by our 
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synthesized beam. In Figure 10 structure sizes are represented by 5 cr contours indicated 
in Figures 3, 5 and 6. We measured the sizes in the north-south and east-west directions, 
and then the values measured in the two orthogonal directions were averaged; a half of this 
average was adopted as angular sizes in Figure 10. The averages should correspond to radii 
if features are spherical. One might think that sizes given by least square fitting by a circle is 
an useful alternative. In fact the circle fitting gives roughly similar results, but occasionally 
produces unreliable values due to irregularity of features; therefore, we do not adopt here 
the circle fitting. The error bars in the vertical axes represent a half of differences between 
measures in the two orthogonal directions, representing degree of deviation from a true cir- 
cle. Observational data points are superimposed on simple models of spherical outflows (see, 
e.g., Fong et al. 2003) with different outflow velocities. Model curves arc fitted to the data 
by a least square; free parameters of the models are only a radius and systemic velocity. At 
a glance, data points of CO J = 1-0 and CO J = 2-1 appear to be fitted well by a model 
with an expansion velocity of 30 km s~^. Although, in the case of CO J = 2-1, data points 
in higher velocity parts (i.e., —2 and 67 km s^^ channels) somewhat deviate from the model 
curve of V^xp = 30 km s"-*^, this is presumably because the features in higher velocity ranges 
are not resolved by our synthesized beam. 

If a spherical outfiow is the case, the broad component region should be identified as an 
AGB wind, because spherical outflows are intrinsically harbored only in late-type stars, but 
not in YSOs. With the expansion velocity of Vexp — 30 km s~^, a dynamical scale of a broad 
component region is estimated to be between 1600 and 3100 yr (on the assumption of distance 
2-4 kpc), which is vahd for a time scale of an AGB circumstellar envelope. However, several 
problems still remains in the spherical outfiow interpretation. First, if the broad components 
really originates in a spherically expanding outfiow, a line profile should be a parabolic or fiat 
top shape (a profile of a optically thin line in a spherically expanding shell is fiat-topped, and 
that of a thick line is parabolic when the telescope beam cannot resolve the envelope, Morris 
& Alcock 1977); however, the line profile of broad components of the ^^CO J = 2-1 line looks 
to be different from both of these profiles, although broad components of other lines seem to 
show parabohc shapes. In fact, the ^^CO J = 2-1 fine profile looks a Gaussian profile rather 
than parabolic or flat top proflles. Second, the angular size of a broad component region 
is larger in the CO J = 2-1 line than in the CO J = 1-0 line; this is very unusual if the 
broad component originates in a usual spherical AGB envelope. To explain this, outermost 
part of the broad component region should be somehow heated to enhance the CO J = 2-1 
intensity. To explain these problems, shock between a central spherical flow and ambient 
material might play a key role. If the shock region have relatively large turbulent velocity, 
the region might produce a high-velocity wing in spectra, possibly exhibiting a Gaussian-like 
profile; this region might also exhibit relatively high temperature due to the shock, causing 
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the larger size of the CO J = 2-1 emission. If this is to the point, the highest velocity part of 
the ^^CO J = 2-1 line should exhibit in a shell-like feature in intensity maps. Unfortunately, 
we cannot confirm such a feature due conceivably to insufficient integration time. As a future 
work, higher sensitivity observations, mapping weak wings, would be important. 

Third, properties of ^^CO J = 2-1 (see bottom panel in Figure 10) are significantly 

different from those of CO J = 1-0 and CO J = 2-1. In fact, the best fit model of a spherical 
outflow to the -"^'^CO J = 2-1 data gives Vc^p ~ 10 km s^^ and Vpcak ~ 36 km s~^; these values 
are largely different from those of CO J = 1-0 and CO J = 2-1. One possible explanation 
for this difference is that a ^'^CO broad component feature is affected by absorption. In fact, 
velocity ranges of strong absorption features seen in ^^CO spectra arc overlapped with those 
of the ^^CO broad component. Although the ^^CO line is relatively optically thin compared 
with the ^^CO line, the ^^CO line might be still affected by absorption if its optical depth is 
very large. Another possible explanation for the ^^CO properties is that a spherical outflow 
has a shell structure expanding spherically, and the inside of the shell exhibits different 
kinematics with the shell (for example, a bipolar flow). Actually, detached shell structure 
has been found in post-AGB stars (see, e.g., Ueta et al. 2005). However, in such a case, 
large number density of molecules, which enable us to detect the ^^CO line, is not expected 
in the hollow of the shell, because post-AGB stars usually do not exhibit significant mass 
loss. This explanation using the hollow structure also conflicts with detection of SiO masers, 
because SiO masers are usually emitted from a near-stellar region. 

Another notable point of a broad component is an intensity ratio of broad component 
emission of the SO and SiO (thermal) lines. According to spectra in Deguchi et al. (2004), 
broad components of the SO and SiO (thermal) lines have about the same strength; however, 
in an 0-rich AGE envelope, the SiO thermal line is expected to be substantially stronger 
than the SO line. Though we did not observed the SiO thermal fine in the present research, 
the SiO line is important as a future work to clarify this problem. At the present stage, 
in our opinion, the spherical outflow interpretation is not conclusive, but also it cannot be 
irrevocably ruled out. 

3.2. Properties of the Narrow Component Region 

A rich set of molecular species of IRAS 19312-1-1950 is one of the central scientiflc 
problems of this source, because the number of detected molecular species is unexpectedly 
large as an 0-rich AGB/post-AGB star, and also because some species detected toward 
the source (i.e., HCN, HNC, CN, CS, IIC3N, and methanol) are rarely found in 0-rich 
AGB/post-AGB stars (see, Deguchi et al. 2004). With results in the present observations. 
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the central star most likely has 0-rich chemistry, for reasons as follows. First, a broad 
component is found in the SO line, but not found in that of the CS line; this indicates that 
SO molecules exist in a central, spherical wind, but CS molecules does not. Second, the broad 
component of SO clearly exhibits an intensity peak at the mapping center (corresponding 
to the intensity peak of near-infrared emission); this suggests that the SO broad component 
originates in the central star. Third, the intensity distribution of the CS emission avoids the 
central region. If we rely on the central 0-rich chemistry, C-bearing molecules, such as HCN 
and CN, are conceivably lying in outer regions, corresponding to a narrow component region. 
In fact, narrow components in molecular rotational spectra of IRAS 19312-1-1950 have raised 
difficulty in evolutional status identification of the source, and have been discussed repeatedly 
in our previous papers (Nakashima et al. 2004; Nakashima & Deguchi 2004; Deguchi et al. 
2004). 

As seen in Figure 7, a narrow component region apparently shows symmetric struc- 
ture elongated in the north-west to south-east direction; this structure appears to exhibit a 
systematic motion. For better understanding of kinematical properties of the narrow com- 
ponent region, we made position-velocity maps in two orthogonal cuts; one corresponds to 
an apparent major axis and another is mutually perpendicular to it (cuts are indicated in 
Figure 7). Position angles of the cuts were determined by eye inspection. In Figure 11 we 
present the position-velocity maps. Systematic variation of structure as a function of radial 
velocity is clearly seen in cuts A; this is particularly prominent in the ^^CO J — 1-0 line, 
exhibiting a feature symmetrically placed with respect to Visr ~ 36.5 km s~^. In the maps of 
CO J — 1-0 and J = 2-1, no features are seen in red-shifted sides (i.e., V\sr < 36.5 km s~^); 
this is due to absorption in corresponding velocity ranges. Features of C^^O J = 1-0 and 
^'^CO J = 1-0 seen in cuts B appear to also exhibit a systematic variation; other lines do not 
exhibits any strong systematic variation in cuts B. In our opinion, the systematic variation 
seen in Figure 11 is explained mainly by a bipolar flow rather than a Kcplcrian rotation disk, 
because the C^^O J = 1-0 and ^^CO J = 1-0 lines exhibits velocity structure variation also 
in cuts B. Though weak rotational motion, coexisting with a bipolar flow, cannot be fully 
ruled out, the rotational motion would not be dominant. Note that Nakashima & Deguchi 
(2004) have also reported a bipolar fiow, detected in the HCO+ J = 1-0 line; they suggested 
that the bipolar axis is lying in the north-east to south-west direction, which is perpendicular 
to the direction of a bipolar flow found in the present observations in the CO lines. We note, 
here, some differences of properties between bipolar flows detected in the HCO^ and CO 
observations. First, a size of the CO bipolar flow is about 3 times larger than that of the 
HCO"*' bipolar flow; this means the HCO"*" bipolar flow is lying relatively closer to the central 
star compared with the CO bipolar flow. Second, an outflow velocity of the CO bipolar flow 
(~ 2 km s~^) is 10 times smaller than that of the HCO+ bipolar flow (~ 20 km s""^); this 
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potentially means that inclination angles of these bipolar flows are different if the flows have 
a same outflow velocity. These differences suggest that kinematical properties of a circum- 
stellar envelope of IRAS 19312+1950 alter with distance from the central star. Otherwise, 
two distinct bipolar flows, which have different scale and different outflow velocities, might 
be individually lying in the envelope (i.e., multiple bipolar flows). 

A relatively large size of a narrow component region compared with a size of a broad 
component region causes a rather difficult problem if we assume that all material has been 
expelled from a single central star. Apparently, a narrow component region is several times 
larger than a broad component region; this potentially means that a dynamical time scale 
of the bipolar ffow is larger than that of the central spherical ffow, because an outffow 
velocity of the bipolar flow is roughly 20 times smaller than that of the central spherical 
flow (of course, this depends on the inclination angle of the bipolar flow). AGB researchers 
widely believe that AGB circumstellar envelopes evolve from a spherical flow to a bipolar 
flow (though some post- AGB stars remain in spherical); the reverse case of this (i.e. from 
bipolar to spherical) has been not expected so far both from theoretical and observational 
view points. One possible solution of this problem is to assume a binary system consisting of 
an AGB and a post-AGB stars, such as OH 231.8+4.2 (rotten egg nebula). OH 231.8+4.2 is 
believed to be a binary system including two evolved stars in different evolutional stages: an 
AGB star (mira type variable with SiO masers) and a post-AGB star (exhibiting a bipolar 
molecular flow) (Sanchez Contreras et al. 2004). In fact, IRAS 19312+1950 is fairy resemble 
with OH 231.8+4.2 in many aspects: a rich set of molecular species, detection of masers and 
a largely extended bipolar flow. Nevertheless, large luminosity expected in a binary system 
causes a problem again: if the binary system is the case, absolute luminosity of the source, 
at least, would be twice of the single-stellar luminosity. A larger luminosity places the object 
further in distance, forcing envelope physical parameters more extreme as an AGB star (see, 
e.g., Deguchi et al. 2004; Nakashima ct al. 2004). Besides, if the binary system includes 
an AGB star, the AGB star would most likely exhibit time variation in infrared luminosity 
(see, e.g., Nakashima et al. 2000); however, recent near- infrared monitoring found that IRAS 
19312+1950 does not show large variation in near-infrared, though small luminosity variation 
cannot be denied (Pujii et al. 2004, in private communications). To pursue this problem, 
more reliable distance is required; SiO masers, detected toward the source, may be helpful 
to determine the distance by measuring the annual parallax if we use a high-sensitive VLBI 
array. 
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4. Summciry 

In this paper, we have reported results of high-angular-resolution radio observations of 
the highly unusual SiO maser source, IRAS 19312+1950, in the CO (J = 1-0 and J = 2-1), 
i^CO (J = 1-0 and J = 2-1), C^^O (J = 1-0), CS (J = 2-1), SO (Jk = 32-2i) and HCO+ 
(J = 3-2) lines with the BIMA array. Main results are as follows: 

1. In line profiles, two kinematical components are seen as reported by previous single-dish 
observations: a broad pedestal component and narrow component. 

2. A broad component region, traced by the ^^CO line, apparently exhibits spherical prop- 
erties with expanding velocity VJgr ~ 30 km s~^; this suggests that a broad component 
originates in an AGB wind. However, several problems still remain; for example, a 
Gaussian-like line profile is seen in the ^^CO J = 2-1 line; the ^^CO J = 2-1 emission 
is more spatially extended than the ^^CO J = 1-0 emission. 

3. A narrow component region traced by CO lines apparently exhibits a bipolar flow, 
which is spatially larger than a broad component region. The dynamical time scale of 
the bipolar flow might larger than that of a possible spherical flow if we observe the 
source nearly in the pole-on direction, causing a conflict with an interpretation by a 
single AGB star. 

4. Intensity distribution of the CS hue avoids the central region of the source, and that 
of a broad component of SO exhibits a small feature, not resolved by our synthesized 
beam, exactly at the mapping center; this means that the central star has 0-rich 
chemistry, and also that C-bearing molecules are conceivably lying at the outer part 
of the nebulosity. 

This research has been supported by the Laboratory for Astronomical Imaging at the 
University of lUinois and NSF grant AST 0228953. This research has made use of the 
SIMBAD and ADS databases. 
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Fig. 1.— Spectra of IRAS 19312+1950 in tlie CO J = 2-1, CO J = 1-0, ^^CO J = 2-1, 
^^CO J = 1-0, C^^O J = 1-0, CS J = 2-1, SO Jk = 32-2i and HCO+ J = 3-2 lines. The 
spectra were averaged over a circle with a diameter of 15"; the averaging circle centered on 
the mapping center. 
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Fig. 2. — Velocity channel maps of the ^^CO J = 1-0 line. The channel velocities are 
indicated in the upper-left corners of each panel. The velocity channels are averaged over 
1 km intervals. The contours start from a 5 cr level, and increase by every 5 a. The 
1 a level corresponds to 4.96 x 10~^ Jy beam~^. The synthesized beam is indicated in the 
lower-right corner. 
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Fig. 3. — Velocity channel maps of the ^^CO J = 2-1 line. The channel velocities are 
indicated in the upper-left corners of each panel. The velocity channels are averaged over 
1 km s~^ intervals. The contours start from a 5 cr level, and increase by every 5 a. The 
1 a level corresponds to 2.24 x 10~^ Jy beam""^. The synthesized beam is indicated in the 
lower-right corner. 
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Fig. 4. — Velocity channel maps of the C^^O J = 1-0 line. The channel velocities are 
indicated in the upper-left corners of each panel. The velocity channels are averaged over 
1 km intervals. The contours start from a 5 o" level, and increase by every 2 a. The 
1 (7 level corresponds to 4.66 x 10~^ Jy bearn"^. The synthesized beam is indicated in the 
lower-right corner. 
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Fig. 5. — Velocity selected channel maps of the ^^CO J = 1-0 line. The channel velocities 
arc indicated in the upper-left corners of each map. The velocity channels are averaged over 
1 km s^^ intervals. The contours start from a 5 cr level, and increase by every 5 cr. The 
1 a level corresponds to 9.93 x 10^^ Jy beam^^. The synthesized beam is indicated in the 
lower-right corners of each panel. The upper, middle and lower panels have different velocity 
steps: 4, 1 and 2 km s~^. 
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Fig. 6. — Velocity selected channel maps of the ^^CO J = 2-1 line. The channel velocities 
are indicated in the upper-left corners of each map. The velocity channels are averaged over 
1 km s~^ intervals. The contours start from a 5 a level, and increase by every 2 a. The 
1 a level corresponds to 4.79 x 10~^ Jy beam~^. The synthesized beam is indicated in the 
lower-right corners of each panel. The upper, middle and lower panels have different velocity 
steps: 4, 1 and 3 km s~^. 
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Fig. 7. — Velocity integrated intensity maps of the CO lines (narrow component). The thick 
and thin contours represent regions in velocity ranges of 35-36 km s~^ (thick) and 37-38 km 
s~^ (thin), respectively. The synthesized beams are indicated in the lower- right corners of 
each panel. The lowest contours correspond to So" level. The steps of the contours and the 
la rms levels of each data are 10 a, 7.02 x 10"^ Jy beam'^ {^^CO J = 1-0), 5 a, 3.39 x 10"^ 
Jy beam-i (^^CO J = 2-1), 5 a, 3.51 x lO'^ Jy beam'^ (^^CO J = 1-0), 10 a, 1.59 x 10"^ 
Jy beam-i (^^CO J = 2-1), and 2 d, 3.30 x IQ-^ Jy beam"! (C^^O J = 1-0). The dotted 
arrows "A" and "B" are cuts used for Figure 12. 
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Fig. 8. — Velocity integrated intensity maps of the CO and SO lines (broad component). 
The left and right columns respectively represent red- and blueshifted wings of the broad 
component. The integration ranges, corresponding to the blue- and redshifted wings, are —6 
- +10, +62 - +79 km s"^ (^^CO J = 1-0), -6 - +10, +62 - +79 km 8"^ (^^CO J = 2-1), 
+21 - +28, +45 - +52 km s'^ (^^CO J = 2-1), and +21 - +28, +45 - +52 km s'^ (SO 
Jk — 32-2i). The synthesized beams are indicated in the lower-right corners of each panel. 
The lowest contours correspond to 5a level. The steps of contours and 1 a rms levels of each 
data are 2 a, 2.41 x 10"^ Jy beam'^ ^^CO J = 1-0), 5 a, 1.16 x 10"^ Jy beam'^ (^^CO 
J = 2-1), 1 a, 8.47 X lO'^ Jy beam"! (^^^O J = 2-1), and 0.5 a, 1.84 x lO'^ Jy beam"! 
(SO Jk = 32-2i). 
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Fig. 9.— Velocity integrated intensity maps of the CS J = 2-1, SO Jk = 32-2i, HCO+ J = 
3-2, and HCO+ J = 1-0 lines (narrow component). The thick and thin contours represent 
maps in velocity ranges of 34-36 km s~^ (thick) and 37-39 km s~^ (thin), respectively. The 
synthesized beams are indicated in the lower-right corners of each panel. The lowest contours 
correspond to 5a level. Steps of contours and 1 a levels of each line data are 2 a, 3.33 x 10^^ 
Jy bcam-i (CS J = 2-1), 2 a, 3.00 x lO^^ Jy bcam-^ (SO Jk = 32-2i), 0.5 a, 4.80 x lO^^ 
Jy bcam-i (HCO+ J = 3-2), and 1 a, 2.20 x lO^^ Jy beam-^ (HCO+ J = 1-0). HCO+ 
J = 1-0 data are taken from Nakashima & Deguchi (2004). The background grey scale 
is a near- infrared composite (J, H, and K bands) image taken by the CI AO camera on 
the Subaru telescope (courtesy of K. Murakawa, M. Tamura, and the CIAO group of the 
National Astronomical Observatory in Japan). 
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Fig. 10. — Size variation of broad component regions of ^^CO J = 1-0 (upper), ^^CO J — 2- 
1 (middle) and ^^CO J — 2-1 (lower) as a function of radial velocity. The curves represent 
models of simple spherical outflows. T4xp represents expanding velocities of the models; the 
unit of the velocity is km s"^. 
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Fig. 11. — Position-velocity maps of CO lines; transitions of lines are indicated in the upper- 
left corners. "A" and "B", indicated in the upper-right corners, mean cuts indicated in 
Figure 7. The cuts A correspond to apparent bipolar axes (see text); the cuts B correspond 
to perpendicular directions to cuts A. Origins in vertical axes correspond to intersecting 
points of cuts A and B (see Figure 7) . The lowest contours correspond to a 5 cr level of rms 
noises. Contour steps are every 7, 3, 1, 5 and 5 (t for ^^CO J = 2-1, ^^CO J = 1-0, C^^O 
J = 1-0, ^^CO J = 1-0 and ^^CO J = 2-1, respectively. 
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Table 1. List of Lines Observed 



moL 


tran. 


rest freq. 


rms 


int. time 


beam 


p. a. 


array 






(GHz) 


(Jy beani"^) 


(hr) 


(arcsec) 


(deg) 






J=l-0 


115.271204 


0.099 


13.9 


4.03x2.76 


3.701 


B, C, D 




J=2-l 


230.538000 


0.479 


12.8 


5.59x2.96 


9.947 


C, D 




J=l-0 


110.201353 


0.050 


12.4 


3.71x2.64 


4.038 


B, C, D 




J=2-l 


220.399000 


0.224 


8.7 


5.33x3.12 


2.265 


C, D 




J=l-0 


109.782160 


0.047 


12.4 


3.84x2.90 


3.956 


B, C, D 


cs 


J=2-l 


97.980953 


0.058 


5.2 


3.51x2.41 


20.29 


B 


HCO+ 


J=3-2 


267.557630 


0.832 


11.5 


7.17x5.88 


-66.05 


D 


SO 


Jk=32-2i 


99.299905 


0.052 


14.0 


3.09x2.64 


41.88 


B 



